Nonextensive theory of dark matter and gas density profiles 
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ABSTRACT 

Pronounced core-halo patterns of dark matter and gas density profiles, observed in relaxed 
galaxies and clusters, were hitherto fitted by empirical power-laws. On the other hand, similar 
features are well known from astrophysical plasma environments, subject to long-range inter- 
actions, modeled in the context of nonextensive entropy generalization. We link nonextensive 
statistics to the problem of density distributions in large-scale structures and provide fundamen- 
tally derived density profiles, representing accurately the characteristics of both, dark matter and 
hot plasma distributions, as observed or generated in simulations. The bifurcation of the density 
distribution into a kinetic dark matter and thermodynamic gas branch turns out as natural con- 
sequence of the theory and is controlled by a single parameter kappa, measuring physically the 
degree of coupling within the system. Consequently, it is proposed to favor nonextensive distri- 
butions, derived from the fundamental physical context of entropy generalization and accounting 
for nonlocality and long-range interactions in gravitationally coupled systems, when modeling 
observed density profiles of astrophysical structures. 

Subject headings: cosmology:theory — dark matter — galaxies: halos — galaxies:structure — plasmas 



1. Introduction 

The analysis of dark matter (DM) and gas 
density distributions in galaxies and clusters is 
presently based on a variety of phenomenological 
models. Following King (1962) relaxed DM halo 
density profiles are well fitted by empirical power 
laws (Burkert 1995; Salucci & Burkert 2000) and 
from N-body simulations by pdm ~ (''/''s)^^(l + 
r/rg)^"^ (NFW), where constitutes a scaling ra- 
dius, chosen to join the asymptotic r-dependence 
(Navarro 1996, 1997). Subsequently a number of 
modifications were proposed (Fukushige & Makino 
1997; Moore et al. 1998, 1999) along with criticism 
of the 'universality' of the NFW-profile (Jing & 
Suto 2000; Borriello & Salucci 2001). Comparing 
the density profiles of DM haloes in cold dark mat- 
ter (CDM) N-body simulations the functional de- 
pendence puM ~ {f/fs)^°'{^ + r/rs)^('^~") (Zhao 
1996) was found to provide good fits to all haloes, 
from dwarf galaxies to clusters at any redshift (Ri- 
cotti 2003), where a is related to the spectral index 
of the initial power spectrum of density perturba- 



tions. Recently a universal density profile for dark 
and luminous matter was suggested (Merritt et al. 
2005). Physically, we regard the DM halo as a 
self-gravitating coUisionless system of weakly in- 
teracting particles in dynamical equilibrium, scale 
invariant from galactic to cluster scales (Burkert 
2000; Firmani et al. 2000; Spergel & Steinhard 
2000). 

On the other hand, the phenomenological 
/3-model pgas /r^)-^/^f^ (Cavaliere 1976), 

where Tc is the core radius, or the double /3- 
model, a convolution of two /3-models with the 
aim of resolving the /3— discrepancy (Bahcall & 
Lubin 1994), provide reasonable representations 
of the hot gas density distribution in galaxies and 
clusters (Xue &Wu 2000; Ota & Mitsuda 2004). 
Physically, /3 corresponds to the ratio of kinetic 
DM and thermal gas energy, assuming values ~ 
2/3. 

Since any astrophysical system is subject to 
long-range gravitational or electromagnetic inter- 
actions, the present situation motivates to intro- 



1 



ducc noncxtcnsivc statistics as theoretical basis for 
both, DM and hot plasma density profiles, uti- 
lized successfully to understand observed core-halo 
structures in astrophysical plasmas (Lcubncr 2004; 
Leubner & Voeroes 2005). In this situation a sin- 
gle parameter k characterizes the degree of nonex- 
tcnsivity or coupling within the system. The cor- 
responding derived power-law distributions consti- 
tute a particular thermodynamic equilibrium state 
(Trcumann 1999), commonly applied in astrophys- 
ical plasma modeling (Leubner & Schupfer 2001; 
Leubner 2002). 

The concept of the Boltzmann-Gibbs-Shannon 
(BGS) thermo-statistics constitutes a powerful 
tool whenever the physical system is extensive, 

1. e. the entropy is additive. This situation holds 
when microscopic interactions and memory are 
short ranged and the environment is an Euclidean 
space-time, a continuous and diffcrentiable mani- 
fold. However, astrophysical systems are generally 
subject to spatial or temporal long-range inter- 
actions making their behavior nonextcnsive. A 
generalization of the BGS entropy for statistical 
equilibrium was introduced from first principles by 
Renyi (1955) and Tsallis (1988), suitably extend- 
ing the standard additivity to nonextensivity. The 
main theorems of the classical Maxwell-Boltzmann 
statistics admit profound generalizations within 
nonextensive statistics leading to a variety of phys- 
ical consequences, sec e.g. Tsallis (1995). Those 
include a reformulation of the classical N-body 
problem (Plastino, Plastino, & Tsallis 1994) or the 
development of nonextensive distributions (Silva, 
Plastino, & Lima 1998; Almeida 2001; Andrade et 
al. 2002) where the duality of nonextensive statis- 
tics, we will focus on, was recognized (Karlin, 
Grmcla & Gorban 2002). Astrophysical applica- 
tions (Plastino & Plastino 1993; Kaniadakis et 
al. 1996; Nakamichi et al. 2002), provided further 
manifestation of nonextensivity in nature. For a 
reformulation in the context of special relativity 
see Kaniadakis (2002). 

2. Theory 

The generalized entropy S{k) characterizing 
systems subject to long-range interactions and 
couplings in nonextensive statistics reads (Tsallis 
1988; Leubner 2004) 



where pi is the probability of the i*^ microstate, 
kB is Boltzmann's constant and the 'entropic in- 
dex' K denotes a coupling parameter quantifying 
the degree of nonextensivity, or cquivalcntly sta- 
tistical correlations, within the system. A crucial 
property of this entropy is the pseudo-additivity 
for given subsystems in the sense of factorizability 
of the microstate probabilities. The transforma- 
tion K = — q) links the k- formalism to the 
Tsallis q-statistics (Leubner 2002). Here k is de- 
fined in the interval — oo < k < cx) where k = oo 
represents the extensive limit of statistical inde- 
pendence and recovers the classical BGS entropy 
as Sb = — fcsX^Pilnpi. Considering two subsys- 
tems A and B the nonextensive characteristics can 
be illuminated in view of the entropy mixing by 
S,{A + B) = S,{A) + S,iB) + S,{A)S,{B)/k, 
a relation consistent with Eq. (1) where the last 
term accounts for the couplings, k < leads to 
an entropy decrease providing a state of higher or- 
der, whereas for «; > the entropy increases and 
the system evolves into disorder. Hence, k can be 
interpreted as a bifurcation parameter measuring 
the two statistical realizations of ordering or disor- 
dering through correlations. A generalization for 
multiple subsystems is discussed by Milovanov & 
Zelenyi (2000). 

Since entropy and probability distributions 
reside physically on the same level the corre- 
sponding generalized energy distributions follow 
as /±(t;) = A± [1 -h v^/iKcr"^)]"" (Silva, Plastino, 
& Lima 1998; Leubner 2004). The superscripts 
± correspond to positive or negative values of k, 
A^ are proper normalization constants, and a de- 
notes the velocity dispersion or mean energy of the 
distribution characterising their width (variance). 
Negative values of k are conveniently introduced 
by changing the sign at appearance, which gener- 
ates a cutoff at V = y/na, see Leubner (2004) for 
details. 

Upon generalization to a spherical symmetric, 
self-gravitating and coUisionless N-body system 
the corresponding steady state phase-space distri- 
bution f{r,v) obeys the Vlasov equation. If the 
particles (stellar system, galaxies) itself provide 
the gravitational potential and f{r,v) is regarded 
as the mass distribution then Poisson's equation 



2 



10" 



„ 10 

T3 



O 

10 



10" 



10" 




10" 



10" 10' 
log(r) [normalized] 



10' 



Fig. 1. — Nonextensive density profiles: for k = 3 

the innermost and outermost curves correspond to 
the /3~ and solutions, respectively. For increas- 
ing K = 4. ..10 both sets of curves converge to the 
central isothermal sphere solution (k = oo) indi- 
cated by dots, p~ from inside and p+ from outside. 



A$ = 47rGp reads 



A$ = A-kG j /(iw^ + ^)d^v, 



(2) 



and represents the fundamental equation gov- 
erning the equilibrium of the system, where 
/(^f^ -|- <I>) = ,f{E) depends on the energy 
only. Commonly, the relative particles energy 
Er = —l/2v^ + ^ is introduced where ^ is the 
relative potential ^ = — $ + $0, chosen to vanish 
at the systems boundary and satisfying Poisson's 
equation as = —iirGp. 

If f{Er) resembles the exponential mass dis- 
tribution function defining the structure of an 
isothermal self-gravitating sphere of gas, in this 
case identical to the phase-space density distribu- 
tion of a coUisionless system of particles. 



f{Er 



Po 



(27ra2)3/2 



exp(- 



) (3) 



then the corresponding density distribution p = 
Po exp(^'/cr^) is found after integrating over all ve- 



locities. Combining with Poisson's equation (2) 
the solution governs the structure of the isother- 
mal self-gravitating sphere (Binney & Tremaine 
1994). The equilibrium distribution Eq. (3) can 
be obtained by extremizing the standard BGS en- 
tropy with regard to conservation of mass and en- 
ergy. 

Since Eq. (3) applies exclusively to a system of 
independent particles we introduce now long-range 
interactions by the generalized entropy functional 
(1). Extremizing Eq. (1) after replacing the en- 
tropy function flnf of an uncorrelated ensemble 
by the generalized functional — k/(1 — f^~^^'^ and 
applying Lagrange multipliers (Plastino & Plas- 
tino 1993) the resulting distribution function reads 
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As previously, the superscripts refer to the 
positive or negative intervals of the entropic in- 
dex K, accounting for less (+) and higher (-) or- 
ganized states and thus reflecting the accompa- 
nying entropy increase or decrease, respectively. 
If we identify f^{Er) again as mass distribu- 
tion, the K-dependent generalized constants 
assure proper normalization and dimension and 
differ for positive and negative definite K-values as 
S+ =Cr(K)/(K3/2r(K-3/2)) and =CT{k + 
5/2)/(/«3/2r(K + l)), where C = po/i2nay/^ and 
r denotes the standard gamma function (Leub- 
ner 2004). The different normalization is caused 
by the interval corresponding to negative k- values, 
which generates an energy cutoff in Eq. (4) leading 
to the constraint — < na^ and restricting 
also the integration limits in velocity space. For 
K ^ oo Eq. (4) approaches the exponential distri- 
bution function (3) defining the density profile of 
the isothermal sphere. 

After incorporating the sign of k into Eq. (4), 
we perform separately for positive and negative 
definite k the integration of (4) over all velocities 
where B^ must be used consistently. The result- 
ing solution provides a modification of the velocity 
space context introduced by Leubner (2004) for 
the density evolution of a system in gravitational 
potential as 



=Po 



1--* 

K a'' 



3/2-K 



(5) 



3 




10 10 10 

log(r) [normalized] 

Fig. 2. — Comparison of the DM nonextensive 
density profile (k = —7, a = 1, solid) with 
the NFW-profile {r, = 0.65, po = 1-2, dashcd- 
dotted). The radial nonextensive gas distribu- 
tion (k ~ 7) is compared with a single /3— model 
(/3 — 0.7, dashed) and the decomposition of a dou- 
ble /3— model {Pcore ~ 0.98, fihaio = 0.55, dotted). 
As illustration a dependence p ~ is added. 



In analogy to the corresponding tandem charac- 
ter in velocity space Eq. (5) generates in a grav- 
itational potential for finite positive values of k 
pronounced density tails, whereas for negative re- 
values the solutions are restricted within the cutoff 
at K = I cP" and k = — oo. 

The duality of equilibria in nonextensive statis- 
tics is manifest in two families, the nonexten- 
sive thermodynamic equilibria and the equilibria 
of kinetic equations, where both are related by 
q =2-q (Karlin, Grmela & Gorban 2002). Since 
q = 1 — 1/k (Leubner 2002) we find for the en- 
tropic index k = —k relating the two families of 
equilibria, where with regard to Eq. (5) k > 
corresponds to the stationary states of thermody- 
namics and K < to kinetic stationary states. The 
limiting BGS state for k = oo is therefore charac- 
terized by self-duality. The nonextensive parame- 
ter K finds also a physical interpretation in terms of 
the heat capacity of a medium (Almeida 2001). A 
system with k > represents an environment with 
finite positive heat capacity and vice versa, for 



K < the heat capacity is negative. Negative heat 
capacity is a typical property of self-gravitating 
systems, see e.g. Firmani et al. (2000). More- 
over, contrary to thermodynamic systems where 
the tendency to dis-organization is accompanied 
by increasing entropy, self-gravitation tends to re- 
sult in higher organized structures of decreased 
entropy. Consequently, the nonextensive bifurca- 
tion of the singular isothermal sphere solution into 
two distributions f^{Er) or p^, respectively, re- 
quires to identify the density profile (5) for pos- 
itive definite k as the proper distribution of the 
thermodynamic state of the gas, whereas the neg- 
ative definite counterpart is associated to the self- 
gravitating DM distribution. For k; — > oo both 
solutions merge at the isothermal sphere distribu- 
tion defined by Eq. (3). 

After elucidating the tandem character of the 
nonextensive context we combine Poisson's equa- 
tion A^' = —A-kGp with Eq. (5) to find a second 
order nonlinear differential equation for for the ra- 
dial density dependence of a spherically symmetric 
gas and DM distribution as 



<fp ^2 dp 



r dr 

AnGn 
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(3/2-n)a2'^ >o 



n'p^dr' 

p2(^)-l/n = 



(6) 



where n = 3/2 — re is introduced, n corresponds 
to the polytropic index n = 1/(7 — 1) of stellar 
dynamical systems where 7 is the adiabatic in- 
dex. Hence, the limit k = 00 is governed consis- 
tently with 7 = 1 by the equation of state of the 
isothermal sphere (Binney & Tremaine 1994) and 
Eq. (4) denotes the nonextensive generalization. 
Multiplying Eq. (7) in the limit k ^ 00 by r^/p 
and combining the first three terms on the left 
hand side as d/d,r{r'^dlnp/dr) recovers the equa- 
tion defining the conventional isothermal sphere. 
Eq. (7) is solved numerically by simulation of 
the corresponding two first order differential equa- 
tions. 

3. Results and discussion 

As natural consequence of nonextensive entropy 
generalization the isothermal sphere profile bifur- 
cates into two distribution families controlled by 
the sign of k and governed by Eq. (7). The self- 
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Fig. 3. — Nonextensive integrated mass of DM 
and gas components (solid) as compared with 
the NFW-profile (dashed-dotted) and the single 
/?— model (dashed) for the gas component. Pa- 
rameters are the same as for Fig. 2. 

gravitating DM component, a lower entropy state 
due to gravitational interaction, resides besides a 
thermodynamic gas component of increased en- 
tropy. Fig. 1 provides the radial density pro- 
file characteristics and k dependences as computed 
from Eq. (7) for both families, DM distributions 
located below and the plasma distributions above 
the limiting exponential solution. With increasing 
entropic index k both braches merge simultane- 
ously in the isothermal sphere density profile, since 
the power-laws (4) converge for k ^ oo to the 
exponential energy distribution (3), correspond- 
ing to n = oo in Eq. (7). Physically finit values 
of K represent long-range interactions and correla- 
tions within the system, whereas the transition to 
«; = oo in the central curve defines the extensive 
limit of statistical independence. Since we focus 
here on the shapes of radial profiles and their phys- 
ical foundations, normalizations arc conveniently 
applied. As r ^ 0, p = po = 1 and the solution 
meets the physically required condition dp/dr = 
in the origin. 

In Fig. 2 we compare one negative (DM) solu- 
tion to Eq. (7) with the NFW model as weU as 
one symmetrically, positive (gas) solution with a 



single /3— model. On small scales the theoretical 
DM density distribution is characterized consis- 
tently with e.g. Firmani et al. (2000) by a shallow 
core of finite density as r ^ 0. Despite a gradu- 
ally tail-like structure of the theoretical profile the 
halo characteristics are similar to the NFW pro- 
file. Changes in the mean energy or variance a 
generates a radial shift of the entire profile, which 
practically corresponds to variations of the scale 
radius r, of the NFW model. As measure of the 
long-range interactions the dimensionless second 
parameter k of the theory controls the shape and 
characteristic mean slope of the profile on inter- 
mediate and large scales, see also Fig. 1. Thus 
both, shallower or steeper slopes (Moore et al. 
1998, 1999), interpreted in the nonextensive con- 
text as consequence of different coupling strengths, 
are accessible. Recently an improved fitting for- 
mula, converging to a finite density in the cen- 
ter, was introduced and compared with the stan- 
dard NFW-profile as well (Navarro 2004). The 
accompanying discussion illustrates critically also 
the potential difficulties arising in simulations of 
the innermost structure of CDM haloes (see also 
Trott & Melatos (2005)). On the other hand, high- 
resolution rotation curve analyses of galaxies are 
consistent only with cored haloes (Blok et al. 2003; 
Gentile et al. 2004) supporting the physical conse- 
quences of the nonextensive statistics formulation 
from observations. The nonextensive gas density 
distribution follows a single /?— model in the core 
but deviates by a halo tail formation. This devia- 
tion can be accurately fitted by a double /3— model 
(a decomposition as in Xuc &Wu (2000) is shown 
by the dotted lines) , indicating that the nonexten- 
sive theory provides naturally a theoretical con- 
text able to solve the /3— discrepancy. 

Fig. 3 presents the radial dependence of the in- 
tegrated mass of DM and gas components for sym- 
metric values of K — ±7 as compared to the NFW 
and /3— model, where on observational grounds a 
20% central gas density fraction po is introduced 
for proper visualization. Consistent with Fig. 
2 the NFW integrated mass exceeds the nonex- 
tensive solution slightly, whereas the integrated 
/3— model is practically identical to the generalized 
entropy approach. 

The dual nature of the nonextensive theory pro- 
vides a solution to the problem of DM and gas 
density distributions of clustered matter from fun- 
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damental physics, where both parameters {k, a) 
admit physical interpretation. The bifurcation of 
the density distribution into a kinetic DM and 
thermodynamic gas branch turns out as natural 
consequence of the theory and is controlled by the 
entropic index k, accounting physically for non- 
locality and long-range interactions in nonextcn- 
sive systems. Due to different correlation prop- 
erties of clustered structures, particular DM and 
gas density profiles might be subject to different 
values of k and a, regulating details of possibly 
non-universal, mass dependent profiles. The the- 
ory was found to reproduce accurately also density 
profiles generated by N-body and hydrodynamic 
simulations, a subject to be addressed elsewhere. 
In conclusion, it is proposed to favor the family of 
nonextensive distributions, derived from the fun- 
damental physical context of entropy generaliza- 
tion, over empirical models when fitting observed 
density profiles of astrophysical structures. 
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